Statistics of soliton-bearing systems with additive noise. 
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We present a consistent method to calculate the probability distribution of soliton parameters in 
systems with additive noise. Even though a weak noise is considered, we are interested in proba- 
bilities of large fluctuations (generally non-Gaussian) which are beyond perturbation theory. Our 
method is a further development of the instanton formalism (method of optimal fluctuation) based 
on a saddle-point approximation in the path integral. We first solve a fundamental problem of soli- 
ton statistics governing by noisy Nonlinear Schrodinger Equation (NSE). We then apply our method 
to optical soliton transmission systems using signal control elements (filters, amplitude and phase 
modulators) . 
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Spatially or temporally coherent nonlinear structures 
(soliton, vortex, breather, domain wall, spiral chemical 
wave, collapsing cavern and many others) play an im- 
portant role in the dynamics and statistics of nonlin- 
ear systems. During the past few decades soliton mod- 
els have arisen in fields as diverse as hydrodynamics, 
plasmas, nonlinear optics, molecular biology, solid state 
physics, field theory, and astrophysics. Presumably the 
most impressive practical implementation of the funda- 
mental soliton concept has been achieved in fiber optics, 
where soliton pulses are used as the information carri- 
ers to transmit digital signal at high bit rates over long 
distances (see e.g. [0J^). In long-haul fiber optic com- 
munication systems, the limitations on the bit rate and 
error-free transmission distance are set mainly by the 
spontaneous emission noise added by in-line optical am- 
plifiers. Existing and future optical transmission systems 
can show no measured errors over long time intervals, 
that makes a direct modeling of the bit error rate (that 
must be less than 10~^) almost impractical. An impor- 
tant role is then played by theoretical methods to evalu- 
ate system performance. Even though the noise is weak 
comparatively to the soliton signal, in general one can- 
not use perturbation approach to obtain the error prob- 
ability because errors occur when soliton amplitude or 
position change substantially due to noise accumulation. 
Though dynamical deterministic properties of many non- 
linear systems have been intensively studied during last 
decades, much less is known about statistics of such sys- 
tems. Typically, consideration of statistical properties 
is limited by assumption of the Gaussian statistics and 
calculations of the variances (note, however, works [^[^] 
where non-Gaussian corrections due to soliton interac- 
tion have been analyzed). Difficulties in studies of non- 
Gaussian statistics in nonlinear systems in particular are 
caused by lack of appropriate mathematical methods. In 



this Letter we present a consistent method to derive the 
probability density function in soliton-bearing systems 
with an additive noise. Our approach here is a modifica- 
tion and further development of a formalism to calculate 
the "probabilities of improbable events" . The method is 
based on finding an optimal fluctuation that provides for 
a maximum of probability under given conditions, techni- 
cally it is a saddle-point approximation in the path inte- 
gral for probability. As a speciflc example, we apply our 
general scheme to the calculation of an error probability 
in fiber-optic soliton transmission. 

We start from the NSE with an additive noise 
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Here ^ is white noise with the correlation function 
{i{ti,xi)£,*{t2,X2)) = D5{ti - t2)5{xi - X2), where D 
is the noise intensity. The equation (^) describes trans- 
mission of the signal along the fiber line, then t is the 
propagation distance and x is time. 

In this Letter, we focus on the problem of a single soli- 
ton distortion by the noise. We assume the ideal soliton 
signal \E'(0, x) = cosh~^ x aXt — Q and examine the prob- 
ability distribution of different distortions of the signal 
at a finite i = T > 0. Another important problem is to 
find the probability to detect "one" at a finite distance 
T provided there was no soliton at t = 0. Solutions 
to these problems are obtained below by analyzing the 
noise-induced fluctuations of ^ around the ideal soliton 
form cosh~^ x. We assume that in the soliton units the 
distance T is large T ^ 1 while the noise power is small 
-D <C 1. More precise condition on D will be formulated 
below. To find the probability to lose digital optical sig- 
nal coded by soliton (elementary "one") t — T, one 
should define a particular measuring procedure that is to 
specify a receiver. For example, the presence of the signal 
a,it = T can be established using the value of the integral 
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J_^dx"^*(T,x)'^{T,x). If the window I is large enough 
then the integral is close to 2 for the soliton cosh~^ x. Er- 
rors are caused by the events with the value of the integral 
essentially smaller than 2. Such rare events are described 
by the tails of the probability density function (PDF). 
The focus of our paper is to develop a regular method to 
calculate such (generally non-Gaussian) tails of the PDF. 
In optical applications there are two leading processes 
which can result in these significant (but rare) deviations 
of the measured energy from its mean value. The first 
process is diminishing of the soliton power, characterized 
by the energy integral Q = da; ^'*(r, x)*(r, x)/2, 
which is equal to unity for the ideal signal. The sec- 
ond process is a shift of the soliton position characterized 
by the integral Y = J^^dxx'^*{T,x)'^(T,x)/{2Q) that 
gives the location of the soliton "center of mass" . For 
the ideal signal, Y = 0. It is clear that when the soliton 
almost leaves the detection window {—1,1}, the integral 
one measures can substantially deviate from 
2. Therefore below we will look for the joint probability 
distribution function 'P{Q,Y). 

We parametrize our signal 5* in the following way: 

^> — rj exp (iPx + ia + ir ) [cosh-\z) + v] , (2) 
z^T]{x~y), dr^Tj'^dt, (3) 

where a, [3, rj, y are soliton parameters that may be ar- 
bitrary functions of time, and we defined the "internal 
time" T. The field v describes the continuous spectrum 
of perturbations on the background of the soliton. An 
important observation is that at T 3> 1 the probabil- 
ities of large deviations are determined by fluctuations 
of the soliton parameters. This is because the discrete 
modes are localized on the soliton and the integral effect 
of a continued in time fluctuation can be significant. On 
the other hand, dynamics of the field v spread its fluc- 
tuations over the whole space. We will show below that 
the influence of the continuous spectrum on statistics of 
the "soft" variables a, [3, rj and y is neglibile in the limit 
DT^ ^ 1. The soliton variables themselves are coupled 
dynamically in a strong way. We first restrict our con- 
sideration by the set a, (3, rj and y, then we will consider 
the continuous spectrum and establish conditions when 
it can be neglected. 

Neglecting the continuous spectrum v we get from (|l|) 



be published elsewhere) that averaging over ^ is equiv- 
alent to the Gaussian average over the new noises with 
zero cross correlations and the following dispersions: 



9t7y = ryC, 9*/? = Ci , 5*2/ - 2/3 = C2 
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plus an equation for the phase dta + ydt(3-\- 0^ — Cs- Here 
the new noises C,{t) . . . (sit) are some spatial integrals of ^. 
These equations has been derived before, albeit without 
a careful definition of the statistics of the noises To 
define the statistics, one needs proper regularization of 
the equations, that can be done by considering a limit of 
a finite correlated case (C(OC(0)) ^ .9(0 with a symmet- 
ric function g. As a result, one may show (the details will 



{at)cm=s{t)D/^, (C3(t)C3(o)) =<5(i)(i2 

(C2(0C2(0)) = 5{t)n'D/12rj', (Ci(i)Ci(O)) = 
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It is interesting that there is a single nonzero average 
(C) = D/2ri which means a systematic increase of the 
soliton amplitude due to the noise: (?y(T)) = 1 + DT/2. 

We are going to study phase-independent quantities 
so a variable can be neglected since it neither enters the 
equations for rj, (3, y nor the noise correlators. We may 
now write 'P{Q,Y) as an average over the noise of the 
solutions of the equations (^ satisfying the boundary 
conditions i]=\,y = [3 = Q&tt = Q and rj = Q, 
y — Y &i t = T. Instead of dealing with the equations, 
it is convenient to take the path integral over arbitrary 
functions 77, (3, y, taking the equations (^) into account 
by corresponding (5-functions which can be rewritten as 
exponentials introducing auxiliar fields /i2- Per- 

forming then averaging over the noise we come to the 
standard Martin-Siggia-Rose formalism: 



dtC{t) 



V{Q,Y) = / VpVTjVyVnVniVfi2exp 



where the effective Lagrangian is as follows 

C = 2i[tidtri/'n - /ii?/2 - nidtP + ii2{dty - 2/3)] 

-D[12^? + {2r^^i,f + {^:^l2/rlf]/{&v) (5) 

Here /i(i),/i2(i) are arbitrary functions on the interval 
(0, T) while Mi(T) = since we do not fix a value of j3 at 
the final moment. From now on we shall be interested in 
values |f? — 1| ^ Dt so we omit the term with the linear 
drift of 77. 

Since we are interested in the events with small prob- 
ability, we calculate the path integral in the saddle-point 
approximation: InV ~ j£saddiodi. The applicability 
condition of the saddle-point approximation is DT <^ 1. 
The extrema conditions that determine the saddle-point 
trajectory (also called instanton or optimal fluctuation) 
can be found from (|^) . Because we are interested in i ^ 1 
one can neglect the field /i2 in comparison with /ii (as 
follows from the relation dtfii ^ ^2) imposing the cor- 
responding condition dty = at < = 0. The resulting 
equations are 

dtv = 2iDn , idtfi = D-q^il/i - Dn^Tj , (6) 
id^y=^DTjfii, a,Vi=0, (7) 

A solution of (|j^) is written via Bessel functions: 

V^{T~t) [Ci Ji/4(k) + C2 J-i/4(«;)]' , (8) 
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where k = A(T — t)^/2 while the relation of the parame- 
ter A to y and Q is found from the boundary conditions. 
The constants Ci and C2 can be expressed via A, for in- 
stance, ^/2C2 = r(3/4)Ai/'*V^. Other fields are easily 
found now from (^[]|). We are interested here in 77 ^ 1 
(though 1 — 77 can be much larger than its rms value), 
then Y ~ XT^. Therefore, considering a region Y ^ T 
we get AT^ <C 1. That procedure corresponds to taking 
only the first terms of the expansion of the Bessel func- 
tions in (j^). Then we get a contribution which is of the 
second order over Y (that is Gaussian as a direct result 
of the applied perturbation procedure with AT^ <C 1): 

in7'(r,g,r)«-— (VQ -1) -RiQ)j^^ (9) 

i?(Q) = 10(1 + 8v/0 + Q)(6 + 3v/Q + (10) 

This joint energy-timing PDF for the NSE is obtained 
for the first time. Remind that it is correct at F ^ T 
and DT"^ < 1 while DT'^ is arbitrary. The PDF is Gaus- 
sian with respect to timing and non-Gaussian with re- 
spect to energy. The most important feature of (^) is the 
consistent analytical confirmation of the empirically well- 
known fact that the dispersion of the timing 4Z?T^/9 is 
much larger than that of the energy {DT) so that the er- 
ror rate of any receiver with an integration window / <C T 
is determined by the timing jitter due to Gordon-Haus 
effect |^,|| . Note that the probability of detecting "one" 
formed from noise (without any soliton initially present) 
is given by a similar instanton solution because an op- 
timal way for a weak noise to create a large signal is to 
grow a soliton. Solving the saddle-point equations with 
the boundary conditions 77(0) = and ri{T) ^ rjf we find 
the probability V^qf) ^ exp{—2r]f / DT) to observe the 
amplitude 77/. Our achievement here is factor 2. 

Next, we examine more sophisticated schemes of opti- 
cal soliton transmission designed to suppress the timing 
jitter. Their general feature is that to compensate for 
the effect of a weak noise it is enough to modifies the 
system only slightly so that an analysis can be done as 
an extension of the above one. 

We consider first the phase modulation which is de- 
scribed by the equation 

-i^t* = a^* + 2|*p*-|-C-e2;2*, (11) 

where the term with e is regarded to be small. It produces 
an additional contribution to the Lagrangian = iiefiiy 
to be added to (||) . Varying the sum C + we get the 
saddle-point equations 

dty = 2P + iD—fi2, dtP = -2ey - —Dtj^ii , 
dtfJ.1 = 2fi2 , dtH2 = -2eAii , (12) 

For r/ and fj, we get the same instantonic equations (|^). 
Below, we will be interested in distances T ^ 1/ ^/e when 
the additional term plays an important role. 



One may show that the evolution of 77 and fi is weakly 
influenced by other degrees of freedom if 7/ <C e^^^^. 
Then we can develop a scheme similar to that in the 
basic case: We examine first the dynamics of 77 and then 
the dynamics of y on its background as a perturbation. 
Note that the physics is different now: The variable y 
oscillates in a parabolic potential while the amplitude of 
the oscillations grows secularly with time. Solving the 
equations we get for T > e^^: 

in7'(T,g,y)«-— (To -1) -i?p(g)-^, (13) 

with RpiQ) = 3(1 + ^/Q + Q)"^ We see that for window 
I > e^-^^^ the fluctuations of both timing and amplitude 
contribute the error rate, that can be calculated using 
the (substantially non-Gaussian) PDF (p^. 

Let us discuss now the role of the continuous spectrum. 
First, it is coupled to the discrete degrees of freedom al- 
ready in the linear approximation because of noise. In- 
deed, if we denote m the fleld conjugated to v (exactly 
in a way fj, fields are conjugated to the discrete variables) 
then the terms proportional to Dmfi appear in the ac- 
tion. A straightforward linear analysis shows that m <ti fJ. 
and those mixed terms can be neglected provided T 3> 1. 
The physical mechanism behind that is the frequency gap 
in the continuous spectrum which makes the mixing non- 
resonant. Second, the continuous spectrum influences the 
soliton parameters via nonlinear interaction. The most 
essential interaction is related to the terms in the La- 
grangian containing /it;^. Note that the fluctuations of 
the continuous spectrum grow with time, this effect is 
related to the noise distributed over the whole space and 
is therefore insensitive to the presence of the soliton. In- 
tegrating over the continuous spectrum, we can find fluc- 
tuation corrections to the reduced Lagrangian associated 
with the nonlinear interaction. A relevant contribution 
to the Lagrangian is ~ {DTYfj? . It has to be compared 
with the bare term D^^. Thus we come to the conclu- 
sion that our scheme is valid if DT^ ^ 1. Note that our 
analysis of the continuous spectrum is not sensitive to 
the presence of the (weak) parabolic potential (basically, 
because of the non-dissipative character of the phase con- 
trol). Therefore the criteria T > 1 and DT'^ < 1 are the 
same for the phase modulation case. 

The most elaborated control scheme that we consider 
in this Letter is when filters and amplitude modulators 
are inserted along the propagation line. This scheme is 
dissipative and it allows one to saturate completely the 
growth of the dispersions of both amplitude rj and timing 
y with an obvious potential for an unlimited propagation 
or information storage ||l|,H,D. We analyze below finite 
fluctuations and discover a new (collapse) mechanism of 
the signal loss which restricts the propagation distance 
or storage time. The propagation equation in this case 
reads (see e.g. Q): 
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where all e's are assumed to be small parameters. The 
ei-term describes an additional amplification necessary 
to compensate the losses due to filtering (e2-term) and 
amplitude modulation (ea-term). Without noise, one has 
a steady soliton with (3 = y = and an amplitude rjs sat- 
isfying €1 — e2?7s/3 + e37r^/1277^. It is linearly stable for 
4e2??s > ^^^sVl^i the condition is assumed to be satisfied 
as well as £3 > 2e2?7f/9 that provides for the stability of 
zero . The expression for the dispersions of the energy 
and timing has been derived before and can be found in 
1^ . Here we describe some properties of the whole joint 
PDF including its time-dependent part responsible for a 
total loss of the signal. The additional terms produce an 
additional effect provided T ^ 1/e, the inequality will 
be implied below. 

An additional contribution to the reduced action 



Ce = iieifi - -ie2firi- - 4ie2/i/3 



(14) 



-e2/ii?? /? - 2ie3/x ( — + 2?; ) —ii3fJ-2y/v 



gives the following saddle-point equations 
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This system is too complicated to solve it analytically yet 
the most important feature can be understood: If y{T) = 
Y is sufRciently large then the amplitude 77 collapses to 
zero in a finite time according to dtf]^ = —e^ir'^/S. For 
V{Q,Y), that means that there is a critical value ^ 1 
so that V{Q,Y) falls into d{Q) if |y| > ¥„■ Of course, 
Yc-[ is a complicated function of Q, ei, £2, £3, that can be 
found only numerically. Below we assume all epsilons to 
be of the same order: £1 ~ £2 '--^ £3 '--^ £. 

Let us examine the region of parameters \Y\ < Kcr, 
Q '--^ 1 that is J/ ~ 1 and 1 — 7/ ^ 1. Then the life time 
of the corresponding instanton can be estimated as £~^. 
Next, we come to estimates 
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So, we can conclude that at |y| < Fcr the stationary part 
of the PDF is as follows 



lnP(Q,r) 



D 



i^^(£2/£l,£3/£l,Q,>^), 
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where -F is a dimensionless function of order unity that 
can be found only numerically by finding the extremum of 
jj" Cdt, which is, evidently, much simpler than a massive 
direct simulations of the noisy NSE. At T > l/£ the PDF 
decays fast in the region \Y\ > ycr, near the boundary 



one can estimate ~\nV{Q,Y) - e{Y - Y„f/D. Thus 
the region Y > practically doesn't contribute to the 
probability of the signal lost. 

The possibility of the collapse leads to the following in- 
teresting and practically important phenomenon. There 
is a finite probability per unit time for the amplitude to 
escape the stability region without returning. This prob- 
ability can be found as a result of the competition of the 
returning terms in the equation for 77 and the noise (^1 
which indirectly influences r] through pumping /3 and y. 
The result is the linearly growing probability of the total 
loss of the signal 
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rexp(-Feoi), F,oi^eyD, 



(17) 



Thus we see that there is a limit time for keeping the 
information. 

The analysis of the continuous spectrum in this case is 
slightly different due to dissipative character of the ad- 
ditional terms. Therefore a saturation of the amplitude 
is observed which of course is £-dependent and tends to 
infinity when £ — )■ 0. Estimation of the continuous spec- 
trum fluctuations give a condition D <C e^^"^ for the above 
scheme to be valid. Note that the conditions of practical 
applicability is more restrictive: (Y'^) ~ D/eie2£3 <C 1 
(otherwise the signal will be lost already at T ~ 1). Note 
that £) <C £'^ is also the applicability condition of the 
saddle-point approximation as is seen from ([T^). 

In conclusion, we have developed a consistent math- 
ematical method to derive the probability density func- 
tions in soliton-bearing systems with additive noise. The 
method is general and powerful enough to enable finding 
probabilities of large deviations in practical propagation 
schemes. 
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